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ABSTRACT 

An XMM-Newton observation of the nearby "pre-cataclysmic" short-period 
{Porb — 3.62 hr) binary QS Vir (EC 13471-1258) revealed regular narrow X-ray 
eclipses when the white dwarf passed behind its M2-4 dwarf companion. The 
X-ray emission provides a clear signature of mass transfer and accretion onto the 
white dwarf. The low-resolution XMM-Newton EPIC spectra are consistent with 
a cooling flow model and indicate an accretion rate of M — 1.7 x lO~ 13 M /yr. At 
48 pc distant, QS Vir is then the second nearest accreting cataclysmic variable 
known, with one of the lowest accretion rates found to date for a non-magnetic 
system. To feed this accretion through a wind would require a wind mass loss 
rate of M ~ 2 x 1O _12 M yr -1 if the accretion efficiency is of the order of 10%. 
Consideration of likely mass loss rates for M dwarfs suggests this is improba- 
bly high and pure wind accretion unlikely. A lack of accretion disk signatures 
also presents some difficulties for direct Roche lobe overflow. We speculate that 
QS Vir is on the verge of Roche lobe overflow, and that the observed mass transfer 
could be supplemented by upward chromospheric flows on the M dwarf, analo- 
gous to spicules and mottles on the Sun, that escape the Roche surface to be 
subsequently swept up into the white dwarf Roche lobe. If so, QS Vir would be 
in a rare evolutionary phase lasting only a million years. The X-ray luminosity 
of the M dwarf estimated during primary eclipse is Lx = 3 x 10 28 erg s -1 , which 
is consistent with that of rapidly rotating "saturated" K and M dwarfs. 



Subject headings: binaries: close — X-rays: binaries — stars: cataclysmic vari- 
ables — accretion, accretion disks 
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Introduction 



Short-period binaries composed of a white dwarf and late-type star are of fundamental 
importance to astrophysics as they are progenitors of cataclysmic variables (CVs) and novae, 
some of which likely evolve to form type la supernovae. Stellar evolution theory predicts that 
these systems are the outcome of a common envelope evolutionary phase in which frictions 
leads to a rapid spiraling down of the orbit of an initial wider binary ( lPaczynskilll976l ) . At the 
conclusion of the common envelope phase, the orbit of the system is expected to decay further 
as a result of angular moment um lost to gravitat i onal radiation and to the magnetized win d 
of the unevolved component ( iKraft et al.l Il962t iPaczyhskil Il967t IVerbunt fc Zwaanl Il98ll ). 
As their separation decreases, the Roche lobe of the secondary is reduced in size and is 
eventually completely filled, initiating the onset of mass transfer onto the white dwarf. 

The timescale for initiation of mass transfer in a post-common envelope binary de- 
pends critically on the angular momentum loss rate. While angular momentum loss through 
gravitational radiation is theoretically well-understood for the stars in CVs, there is no com- 
prehensive theory of spin-down through magnetized winds that is expected to dominate the 
evolution of systems with periods above the CV period gap of > 3hr. It depends on the wind 
mass loss rate and t he large scale stellar magnetic field (e.g. IWeber fc David Il967t iMestel 
19681 ; iKawalerl Il988h . Both are extremely difficult to measure for late- type stars and are 
especially uncertain at the very rapid rotation rates of close binaries where magnetic proxies 
such as X-ray emissi on show saturation and supersatura tion effects that appear to affect 



spin-down rates (e.g. Ivanova &: Taam 2003 ; Barnes! I2003I ). 



In this paper, we present XMM-Newton observations of the intriguing DA+M2-4 3.6hr 
eclipsing binary QS Vir (formerly more commonly known as EC 13471-1258) that was first 
discovered as a whit e dwarf-M dwarf binary in the Edinburgh-Cape blue object survey 
( Kilkenny et al.lll997| ). QS Vir has gained considerable recent attention as a valuable diag- 
nostic of close binary evolution and angular momentum loss. Based on an extensive analysis 
of high-speed and multi-colour photom etry, to gether with UV HS T STIS and visible light 
spectroscopy, lO'Donoghue et al.l ( 120031 . see also lKawka et al.ll2002l ) concluded that the sys- 
tem filled its Roche lobe and was probably undergoing very weak mass transfer and suggested 
the system might be a hibernating cataclysmi c variable. It was detected in the ROSAT All 



Sky Survey (RASS) at 0.14 ± 0.02 count s 1 (jVoges et al.lll999l ) implying an inconspicuous 



X-ray lumi n osity Lx = 5 x 10 erg s — much too low for an accreting CV. More recently, 



Qian et al.l ( 120101 ) inferred the presence of a giant planet with ma s s QAM .j Mpit e.r in a 7.86 yr 
orbit based on long term period variations, although iParsons et al.l (120101 ) ruled this out from 
more complete eclipse monitoring. Instead they find a third body in a highly elliptical orbit 
best fits the more complete data, noting that this explanation for the secular orbit change 
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is also not without problems. iRibeiro et al.l (120 lOf ) found evidence for materi al possibly as- 



sociat ed with wind accretion within the Roche lobe of the white dwarf, while iParsons et al. 



( 120 111 ) interpreted phase-specific absorption signatures as a stellar prominence originating 
from the M dwarf. The latter also found the white dwarf to be slowly rotating, indicating a 
pre-cataclysmic evolutionary status. 

Here, we report the discovery of deep, narrow X-ray eclipses locating the dominant 
source of X-rays on the white dwarf rather than the corona of the late- type star and implying 
ongoing accretion. QS Vir has one of the lowest accretion rates of kn own non-magnetic 



syste ms. This, together with its eclipsing nature and proximity (48 pc; lO'Donoghue et al. 



20031 ). renders QS Vir an interesting object for the study of CV evolution. 



2. Observations and Data Analysis 



QS Vir is an eclipsing binary with an inclination of 74 ± 2° and an orbital period of 
Porb = 0.15074 ±0.00004 days. The component spectral types are classified as DA and M2-4 
with masses M wd = 0.78 ± 0.04 M p, R md = 0.011 ± 0-01 ^ and M md = 0.43 ± 0.04 M ffi 
R md = 0.42±0.02 Rr^ respectively fcawka et ali booi lO'Donoghue et alil2003l : Irlowell et al" 



20101 ; IRibeiro et al.ll2010l : iQian et al.ll2010f ). 



QS Vir was observed usi ng the instruments on b oard the XMM-Newton satellite (IJansen et al 



200l|). Both RGS detectors 



the EPIC-pn (IStriider et al. 



den Herder et al.ll200ll ) operated in "Spect roscopy" mode, while 
200ll) and bothMOS jTurner et al.ll200l[^ cameras operat ed in 



"Small window" mode using a medium filter. The OM telescope (IMason et al.ll200ll ) was 
used in "Fast Mode" with the U filter, but a hardware glitch placed the source outside of 
the photometry window and severely compromised the data. A summary of the observations 
is reported in Table [TJ Here we concentrate on the EPIC observations; RGS spectra were 
of low signal-to-noise ratio, although a number of weak emission lines attributed to O VIII, 
Fe XVII, Ne IX an Ne X could be discerned. 



2.1. X-ray and Optical Photometry 

Data were reduced with the XMM-Newton Science Analysis System software version 
7.0.0, updated with the latest calibration files. The EPIC-pn background-subtracted light 
curve with a 100s bin size is showed in Fig. [TJ Apart from the sizable flare 20 ks into the 
observation, which is likely due to coronal activity on the red dwarf, the most conspicuous 
feature of the light curve is the presence of deep eclipses of duration ~ 800s and a period 
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qf abo ut 13 ks. This eclipse period and duration matches that found by lO'Donoghue et al. 
( 120031 ) in their UV-optical study of QS Vir, and we attribute the eclipses to the UV-bright 
white dwarf passing behind the red dwarf. At mid-eclipse, the X-ray emission declines to 
about 15% of its uneclipsed intensity, indicating that ~ 85% of the X-ray signal originates 
from the vicinity of the white dwarf and dominates the quiescent coronal X-ray intensity of 
the M dwarf. 

OM frames were analyzed using custom software designed to extract the residual spill- 
over signal in the photometry window from point spread function wings. While not photo- 
metric, sufficient signal was recovered to verify that the X-ray and [/-band eclipses coincided 
with each other. 

The EPIC-pn light curve phase- fold ed on the ephemeri s for o ptical mid-eclipse derived 
from observations between 1992-2002 by lO'Donoghue et al.l (120031 ) succeeded in co-aligning 
eclipses in phase space. However, the X-ray eclipse preceded the optically-derived one by 
about 120s. T his is a much larger offset than expected based on the period drift apparent 
in Figure 3c of lO'Donoghue et al.l ( 120031 ). indicating that the orbital period had shortened 
somewhat in the interval between the ephemeris epoch and our 2006 observation s. This sec- 
ular p eriod shortening wa s confi rmed recently by the later optical observations of lOian et al. 
( 120101 ) and iParsons et al.l ( 120101 ) . The ephemerides of these works are similar for the epoch 
of our observations and adopting arbitrarily that of the former we find our X-ray and [/-band 
eclipses are centered on the expected primary eclipse. The phase-folded X-ray light curve 
based on this ephemeris is illustrated in Fig. [2J 

Each cycle of the phase-folded X-ray light curve was fitted with a trapezoidal profile to 
locate the times of ingress and egress. No significant variability in the eclipse mid-time or 
duration are present, and the time of second to third contact is r 2 3 = 795 ±40s, in reasonable 
agreem ent with the value of 835 ± Is obtained from optical photometry by lO'Donoghue et al. 
( 120031 ). This indicates that the eclipsed X-ray source is close to or coincident with the white 
dwarf photosphere. 



2.2. EPIC Spectroscopy 

EPIC pn source photons were extracted from a circular region with a 35" radius, while 
background was estimated from a surrounding annulus with an area of approximately ten 
times the source region. Intervals showing particularly high background were excluded from 
the analysis. The extracted EPIC pn spectrum, rebinned such that each channel contains 
at least 35 counts, is illustrated in Figure [3J 
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Mukai et al.l (120031 ) and iPandel et al.l ( 120051 ) have found that the X-ray spectra of CVs 
fall into two broad categories: optically-thin spectra resembling coolin g flows; and spectra 
with a hard continuum exhibiting signatures of a photoionized plasma. IMukai et al.l (120031 ) 
noted that the difference is likely to be associated with the accretion rate per unit area, with 
photoionization spectra belonging to magnetic systems with channeled accretion streams 
covering a small surface area, and optically-thin spectra coming from the accretion boundary 
layer of non-magnetic systems. 

Cursory inspection of the EPIC spectrum of QS Vir reveals an optically-thin form not 
unlike that seen from active stellar coronae. This is reinforced by the few emission lines that 
could be discerned from the RGS spect ra. We performed param eterized spectral modeling 



within the XSPEC program version 12 (IDorman fc Arnaud 



20011 ). using multi-temperature 



optically-thin plasma radiative loss models and a cooling flow model. A fixed absorption of 
Nh = 10 18 cm -2 , corresponding to negligible interstellar X-ray absorption for this nearby 
object, was adopted for both models. For the multi-temp erature plasma, w e e mployed a two 



temp erature MEKAL model (a melange of the models of lMewe et al.lll985l and lLiedahl et al. 



19951 ) in which the temperatures were allowed to vary freely but the metallicity parameter was 
forced to assume the same value for both components. The best-fit model had a minimum 
value of the reduced \ 2 statistic of Xo = 1-42, with 520 degrees of freedom (dof). For the 
cooling flow, we employed the isobaric flow model with emi ssion measure verses temperature 
scaled inversely to the cooling time at that temperature (IMushotzky fc Szymkowiaklll988l ) 
implemented as the MEKAL-based MKCFLOW model. No redshift was included, and Xo = 
1.23 with 520 dof was obtained with metallicity, minimum and maximum temperatures, 
and normalization as free parameters. The parameter estimation results are reported in 
Table [2] and the best-fit model is illustrated in Figure |3j At face value, the cooling flow 
model describes the data slightly better than the two-temperature one, though we caution 
that systematic uncertainties in the model and instrument response are not included in the 
analysis. 



3. Discussion 
3.1. Accretion 

The detection of clear, sharp X-ray eclipses phased with the optical eclipses provides an 
unambiguous indication of ongoing accretion in the QS Vir system. From the EPIC spectrum, 
we find the total X-ray luminosity in the 0.5-10 keV range to be Lx = 2.1 x 10 29 erg s _1 . 
T he source exhibits n o n-flari ng variability at a level of about 20%. Based on the parameters 



of lO'Donoghue et al.l (120031 ). the bolometric luminosity of the M dwarf is L^i = 5.6 x 
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10 erg s , and the ratio of X-ray to bolometric luminosity is Lx/Lboi = 5.4 x 10 — 
comparable wit h the canonical emp irical "saturation" value of 10~ 3 found for the most 



active stars (e.g. iPizzolato et al.ll2003f ) and typical of the pre- cataclysmic bin aries and rapidly 



rotating M dwarfs in the solar neighborhood studied by iBriggs et al.l ( 120071 ). 



pi 
E 



If QS Vir is a semi-detached system, the accretion rate obtained from the cooling flow 
model, M = (1.69 ± .05) x 1O _13 M yr _1 , is among the lowest ever found for a CV. 



Lynden-Bell fc Pringld (119741 ) showed that approximately half of the gravitational energy 
of the accreting gas is liberated through optical and UV radiation in a viscously-heated 
accretion disk and the other half is dissipated largely in the form of X-rays in a "boundary 
layer", where the disk material is decelerated from its Keplerian velocity to the rotation 
velocity of the white dwarf. The expected X-ray luminosity from a disk boundary layer is 



then L BL ~ GM wd M /2R wd . Using R wd = O.Olli?© and M wd = 0.78M Q (lO'Donoghue et al. 



20031 ). the boundary layer luminosity is Lbl = 7 x 10 erg s , which is reasonably consistent 



with that observed. If the radiative dissipation of viscous heating amounts to a similar 
luminosity in the UV and optical, it would be inconspicuous in comparison to the wh ite 
dwarf luminosity of L wd = 1.7 x 10 31 erg s" 1 (T e// = 14220 K; lO'Donoghue et al.l l2003h as 
expected. 



3.2. Roche lobe overflow or just a wind? 

3.2.1. Difficulties with Roche lobe overflow 

Accreting gas could originate eithe r from the wind of t he re d dwarf, or from Ro che 
Lobe overflow. lO'Donoghue et al.l (120031 ) . iRibeiro et al.l (|2010[ ). and IPar sons et al.l (120111) all 



found evidence for material within the Roche lobe of the white dwarf. 



O'Donoghue et al. 



( 2003 ) suggest this might be associated with an accretion stream, while Ribeiro et al.l ( 201oh 



favored wind accretion because p hotometric modeling suggested the M dwarf did not fill its 
Roche Lobe. iParsons et al.l (1201 if ) suggested the presence of significant prominence material. 



An additional problem with Roche lobe overflow for QS Vir is that this should lead to 
formation of an accretion disk. At low accretion rates a disk can be more optically-thin and 
still render relatively strong disk emission lines, but no observable signatures of a disk have 
been detected. Since our observed mass accretion rate exceeds the total solar wind mass 
loss rate (~2x 1O~ 14 M ) by an order of magnitude, one crucial question to decide between 
accretion mechanisms is whether the mass loss rate of the M dwarf component of QS Vir is 
sufficient to supply this flow through a wind alone. 
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3.2.2. Insufficient wind supply's 



QS Vir is in some respects similar to the "pre-polars" — magnetic (B > 10 6 G) pre-CVs 
that do not appear to fill their Roche lobes yet show signs of accretion activity at low mass 
transfer rates. Accretion on these systems is thought to be through magnetic capture of the 



secondary wind with a wi nd accretion e fficiency of the order of 100% (e.g.. ISchwope et al. 



20021 ; ISchmidt et aDboQsl : see also, e.g., iLi et al.lll99l IWebbink fc Wickramasinghell2005h . 
Accretion rates in the range M = 5 x 10~ 14 -3 x 10 _13 M Q yr _1 with no apparent depen- 
dence on orbital per i od or secondary spec t ral type have been foun d for several systems (e.g. 



Schwope et al.ll200a ISchmidt et al.ll2005l 120071 Vogel et al.ll201ll ). These rates are signifi 



cantly higher than wind accretion on non-magnetic pre-CV systems where magnetic capture 
is not expected to operate. Total mass loss rates for M dwarfs in close non-contact binaries 
with non-magnetic white dwarfs of ~ 10~ 14 -10 -15 M Q yr - 1 have been estimated based on th e 



gravitational settling time of wind-accreted metals (e.g. iDebesI 120061 ; iTappert et al.ll201ll ). 
Schmidt et al.l ( 120071 ) argue that the pre-polar accretion rates could provide the first realistic 
measurements of stellar mass loss at the cool end of the main sequence. Si ngle M dwarf mass 
l oss ra tes have otherwise been notoriously difficult to pin down; see, e.g., IWargelin fc Drake 
( 120021 ) for a detailed discussion. IWargelin fc Drake! ( 120021 ) used astrospheric change exchange 
X-ray emission to place an up per limit of M < 3 x 1O _13 M yr _1 for the moderately ac- 
tive M5.5 dwarf Proxima, while IWood et al.l (120021 ) found M < 4 x 1O~ 15 M yr -1 from Lya 
observations. Both of these are reasonably consistent with the rates estimated for pre-polars. 

The wind mass loss rate required to feed the observed 2 x 10~ 13 M Q yr -1 accretion rate 
on QS Vir depends on its accretion efficiency. In absence of detailed magnetohydrodynamic 
accretion simulations this efficiency m ust be considered un known. Exploratory MHD wind 
models for detached pre-CV systems by lCohen et al.l ( 120121 ) find accretion efficiencies depend 
quite sensitively on the WD and donor magnetic field strengths and on their orbital separa- 
tion. We can rule out QS Vir being a pre-polar with a strongly magnetized WD based on 
the absence of any Zeeman broadening or splitting of the Balmer l ines or cyclotr o n em ission 
features in th e spec tra illustrated by lO'Donoghue et al.l (120031 ). iRibeiro et al.l ( 120101 ) and 
Parsons et al.l (120 111 ). The cores of the Balmer lines in QS Vir are dominated by emission, 
which would obscure any small Zeeman splittings and defines an upper limit to the possible 
field streng th. Based on the ±5 A width of the emission cores in the spectra illustrated in 
Figure 1 of iParsons et al.l ( 1201 ll ). and using Equati on 11 from IWickramasinghe &: Ferrario 
( I2000l ) for a white dwarf temperature of 14200 K ( lO'Donoghue et al.ll2003l ). we estimate 
B < 700kG. Adding some allowance for Stark broadening hiding a slightly higher field, we 
find a conservative limit of B < IMG. Such a limit is also consistent with the absence of 
cyclotron emission feat ures, which are prom i nent in pre-polars with even lower accretion 
rates than QS Vir (e.g. ISchmidt et al.l 120071 ). ICohen et al.l (120121 ) did find one pre-CV con- 
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figuration for a WD field B = 10 5 G in which the accretion rate approached 100% of that of 
a single M star wind by an efficient "siphoning" of plasma from the hemisphere facing the 
WD, but this finding needs to be confirmed by more detailed simulations. The field would 
need to be of order 10 7 G to capture all the wind of the M dwarf, and the accretion efficiency 
is then unlikely to be close to 100% for QS Vir. 

An order of magnitude es t imate of the pure hydrodynamic accretion efficiency can be 
made from the iBondi fc Hoyld (119441 ) recipe 



A7rG 2 M wd 2 p(r) 



(1) 



where the wind density, p(r), as a function of radial distance r in the vicinity of the white 
dwarf at orbital separation, a, for a wind mass loss rate M at velocity v is 



M 



P 



AnvR 



(2) 



rd 



Combining Eqns. [T] and El and assuming the wind velocity is approximately the escape 
velocity, v ~ y^2GM rc i/R r d, (as is the case for the solar wind), the accretion rate is 



RrdMu 

2aM, 



rd 



M 



(3) 



where the first term on the right is the accretion efficiency. The orbital separation, a = 
[P% rb G(M W d + Mrdj/^w 2 ] 1 / 3 , from the parameters in §2] is 1.26i? and the accretion efficiency 
is 0.09. Considering the accretion rate derived for QS Vir from the one-dimensional cooling 
flow model is likely a lower limit, the wind mass loss rate required to feed the accretion is 
then at least M~2x 1O _12 M yr _1 — an order of magnitude larger than implied for the M 
dwarfs of the pre-polars. For a typical velocity of 600 km s -1 , this wind would require ~ 1% 
of the stellar bolometric luminosity to drive which we consider uncomfortably high. 



Qian et al.l (120101 ) interpreted a residual period decrease of QS Vir with respect to their 



ephemeris solution that included a giant planet as angular momentum loss from a wind with 
an even higher mass l oss rate of 7.2 x 10~ n M(n yr -1 . Massive winds, up to 10 _10 M Q yr _1 , 
were argued earlier by iMullan et al.l (119921 ) based on ther mal escape arguments togeth er with 
weak evi dence of infrar e d win d free- free emission (see also lBadalyan fc Livshitslll992l ). More 
recently, IVidotto et al.l (120111 ) predicted a similarly massive wind for the M4 dwarf V374 Peg 
from magnetohydrodynamic modeling, with a ram pressure 5 orders of magnitude larger 
than the solar wind. However, such massive winds can be readily dismissed on energetic 
gro unds. For a termi nal velocity of ~ 500 km s _1 , the kinetic power of the wind suggested 
by iQian et al.l (120101 ) is 6 x 10 30 erg s -1 — an implausibly large fraction (10%) of the M 
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dwarf luminosity. Similarly, the wind of IVidotto et al.l (1201 1 ) corresponds to an unphysical 



10 erg s of kinetic energy wh en compared with the underlying stella r bolometric 



luminosity of L bo i = 3.6 x 10 31 erg s _1 fjMorin et al.ll2008l ). iLim fc White! (119961 ) had pointed 
out earlier that such massive winds would be opaque to radio emission which is commonly 
observed from active late-type stars, and deduced that mass loss rates can be no larger than 
1-2 orders of magnitude above the solar rate. For QS Vi r, the need for such a massive wind 
is obviated by the additional period measurements of iParsons et al.l (120101 ) . who instead 
invoked a third body with a minimum mass of 0.05M Q i n an ellipti c al orb it, together with 
the activity-relate d quadrupole couplin g mechanism of lApplegatd (119921 ). to explain the 
period variations. IParsons et al.l (120101 ) noted evolutionary difficulties for the presence of 
such a body, though it and the secular period evolution are not of any obvious relevance for 
understanding the pre-CV evolutionary status, parameters and accretion rate of QS Vir. A 
further unexplored connection cannot of course be ruled out. 



The kinetic power of the solar wind, 



10 27 erg s \ is similar t o the mean observed 



radiative power of the corona (e.g. iPeres et al.l 120001 ). ISchwadron et al.l (120061 ) argued that 
the energy available to drive the solar wind scales with the basal open magnetic field strength. 
Combining this with an observed, almost linear, correlation between X-ray radiance and 
surface magnetic flux ( Pevtsov et al. 20031). the expe c tation would be for stellar wind mass 
loss rate to scale with X-ray luminosity. IWood et al.l (120021 ) estimated mass loss rates for a 
handful of different stars of different activity using astrospheric Lya absorption, and found 
M oc p 1 ^.- 15=l=0 - 20 ) where Fx is the surface X-ray flux. Scaling to the observed surface X- 
ray flux of QS Vir yields M ~ 5 x 10~ 12 M Q yr _1 — consistent with our estimated accretion 
requirement. Howe yer, there also exist no measurements for stars at the saturated activity 
level of QS Vir and IWood et al.l (120021 ) caution against extrapolating the relation. Proxima 
also does not support the relation, with a wind upper limit from Lya absorption lying an 
order of magnitude below the predicted value. 

In summary, for an assumed wind accretion efficiency of the order of 10%, the prepon- 
derance of more direct assessments of the likely wind mass loss rates of active M dwarfs 
like the secondary of QS Vir, in particular those for the pre-polars, suggest it would not be 
sufficient to feed the observed accretion rate. 



3.2.3. Chromospheric mass supply in a rare evolutionary phase? 



An alternative to catastrophically energetic wind or accretion e fficiency requir e ments 
is surface activity that does not lead to mass loss in single stars. IParsons et al.l (j201l[ ) 
suggested prominence activity could be responsible for material within the Roche lobe of 
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the white dwarf. Such activity is stochastic and might be expected to result in highly 
variable X-ray emission. The X-ray luminosity found here from XMM- Newton observations, 
(2.1 ± 0.4) x 10 29 erg s _1 , is about a factor of two lower than that measured in the RASS 
((4.9 ± 0.7) x 10 29 erg s~ 4 ) — perhaps suggestive of moderate accretion variability though 
flaring on the M dwarf cannot be ruled out. 

A more constant, lower lying, source is upward chromospheric spicule and mottle flow. In 
the solar chromosphere and transition region this flow co rresponds to a mass flux two orders of 
magn itude larger than is launched in the solar wind (e.g. ICamposlll984l ; iTsiropoula fc Tziotziou 
20041 ) . The flow has a typical scale height of ~ 5000 km and falls back to the solar surface, 
but might more easily succeed in escaping the Roche surface of a star on the verge of Roche 
lobe overflow to become readily unbound to the M dwarf and available to be swept up by 
the lobe of the white dwarf. Such a mass flow source w ould only become available for a near 
lobe-filling system. Following, e.g.. iKing fc Kolbl (119951 ). the timescale for evolution through 
such a chromospherically-fed phase is 



h r h J 



Tch 



(4) 



Rrd2J 

where h c h = kT c hR rd 2 j [iGM rd is the chromospheric scale height, J the system angular mo- 
mentum and J its time derivative. For a chromospheric temperature of 10 4 K, the ratio 
hrh I ' Rrd = 4.3 x 10~ 4 , thou gh for spicules this height i s an o rder of magnitude greater. Adopt- 



in g the ISkumanichl ( I1972h- based IVerbunt fc Zwaanl (119811 ) spin-down prescription modified 
by Rappaport et al. (jl983 ) 



3 = -3.8 x lO 3O M rdJ R 4 



R 



rd 



Rp) 



with 7 = 2 and where 



J 



2tt \ a 2 M rd M wd 



2tt 



P 



orb 



(5) 



(6) 



, Porb J M rd + M w d 

we find a timescale t c k — 5 x 10 4 yr for the chromospheric scale height and 5 x 10 5 yr for 
spicules. The implication of this short timescale is that chromospheric mass supply would 
occur only in a rare, fleeting episode of evolution. 



4. Conclusions and the nature of QS Vir 



X-ray spectrophotometry of QS Vir betrays on going accretion at a ra te of 2 x 10 _13 M Q yr~ 
The slow rotation of the white dwarf diagnosed by iParsons et al.l (120111 ) argues against a hi- 
bernating system, and QS Vir is most likely at the first onset of significant accretion. QS Vir 
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then represents a valuable object to probe the very earliest stages of CV evolution. At face 
value, the observed accretion rate implies Roche lobe overflow and a semi-detached system. A 
pure wind accretion alternative would imply a wind mass loss rate of M ~ 2 x 10~ 12 M Q yr _1 
for a Bondi-Hoyle accretion efficiency of 10%. For a typical velocity of 600 km s -1 , this 
wind would require 1% of the stellar bolometric luminosity to drive. While plausible, we 
consid er such a mass i ve an d energetic wind unlikely. If QS Vir is still detached, as the anal- 



ysis of iRibeiro et al. (120101 ) implies, accreted mass could be supplied by prominence activity 
(iParsons et al.ll201ll ). though such stochastic activity would likely lead to highly variable 
X-ray emission. The X-ray luminosity observed by XMM-Newton is a factor of two lower 
than that seen in the RASS and weakly supports this idea. We further speculate here that 
upward chromospheric mass flow analogous to solar spicules and mottles that escapes the 
Roche surface to be swept up by the white dwarf would present a more constant mass sup- 
ply. Such a supply would become available for about million years or so prior to Roche lobe 
overflow and QS Vir would then be in a very interesting, fleeting evolutionary phase. 



We are indebted Boris Gansicke for very useful comments and discussion of low accretion 
rate CV-like systems. MM was supported by NASA Grant NNG06GD34G during the course 
of this work. JJD and VK were funded by NASA contract NAS8-39073 to the Chandra X-ray 
Center (CXC) and thank the CXC director, Harvey Tananbaum, and the CXC science team 
for advice and support. 
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Table 1. Summary of the observations 



Instrument 


Exp. Time 
[ksec] 


Mode 


Filter 


Obs. Start Time 


RGS1 


109.5 


Spectroscopy 




2006-01-25@18:13:37 


RGS2 


109.5 


Spectroscopy 




2006-01-25018:13:42 


M0S1 


109.0 


Small Window 


medium 


2006-01-25@18:14:51 


MOS2 


108.9 


Small Window 


medium 


2006-01-25@18:14:51 


pn 


109.3 


Small Window 


medium 


2006-01-25@18:20:24 


OM 


17.6 


Fast Mode 


U 


2006-01-25@18:20:22 


11 


17.6 


Fast Mode 


U 


2006-01-26@00:53:50 


ii 


21.8 


Fast Mode 


U 


2006-01-26007:57:18 


ii 


21.7 


Fast Mode 


u 


2006-01-26@14:28:41 


ii 


21.7 


Fast Mode 


u 


2006-01-26@20:57:59 
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Table 2. Summary of the parameter estimation results 



COOLING FLOW 


2-TEMPERATURE 


Parameter 


Best-fit 


Parameter 


Best-fit 


N H (cm- 2 ) a 


10 18 


(cm- 2 ) a 


10 18 


low T (keV) 


0.08 ±0.15 


kTi (keV) 


0.67 ±0.01 


high T (keV) 


7.24 ±0.17 


kT 2 (keV) 


3.58 ±0.08 


Metallicity b 


0.44 ±0.05 


Metallicity b 


0.48 ±0.03 


M (lO" 13 M /yr) 


1.694 ±0.046 


EMi (10 52 cm~ 3 ) c 


0.26 ±0.03 






EM 2 (10 52 cm~ 3 ) c 


2.08 ±0.04 



a Parameter held fixed. 

b Expressed as a f ractio n of the solar metal abundances listed by 
Anders fc Grevessd (Il989l ) 



= Assuming a distance of 48 pc ( lO'Donoghue et al.l 120031 ) 




Fig. 1. — EPIC-pn background-subtracted light curve illustrated with a 100 sec bin size 
a function of time in seconds since the beginning of the observation. 
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t 1 1 1 1 1 1 1 1 1 1 1 r 




-0.4 -0.2 0.0 0.2 0.4 

Phase 



Fig. 2. — EP IC-pn back g round subtracted lightcurve phase-folded according to the optical 
ephemeris of iQian et al.l ( 120101 ) . Each phase is offset upward with respect to the previous 
one by 0.5 count s" 1 . 
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Fig. 3. — EPIC-pn spectrum together with the best-fit MKCFLOW cooling flow model. 



